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RESEARCH SUMMARY 


Predicting fire behavior in nonuniform fuel arrays is a problem 
requiring: 


ils A method of assessing fuel nonuniformity, 
Die A method of simulating fuel nonuniformity, and 
S. An algorithm governing fire spread through a simulated array. 


Satisfying these requirements is the objective of this paper. The main 
concept is built around partitioning the fuel into a honeycomb array. 
Each cell is described independently according to its bulk fuel para- 
meters (depth, load, average particle size, etc.). Field assessment is 
designed to meet the requirements of simulation. An algorithm simulates 
fire spread through the array by coupling predictions of heat flowing 
from a burning cell to predictions of the heat required for ignition of 
the adjacent cells. Ignition is allocated to the cell offering the least 
requirement for heat. Consequently, the fire moves nonuniformly through 
the array taking advantage of the path of least resistance. Methodology 
is emphasized. 


A simulated fire is initiated from a line source. Distortions in 
the propagating front result from fuel nonuniformities giving rise to a 
distribution of rates of spread rather than a single value. Analysis 
is appropriate for an assessment of a distribution of the fireline 
intensities. 


Examples are given for slash, residue after tree harvest, and a 
mixture of grass and sagebrush. Nominal windspeeds of 0 and 2 mi/h 
were chosen for the purpose of illustrating the technique for handling 
nonuniformity. 


Comparisons show that the previous alternative of combining all 
fuel to an average depth and load does not allow the land manager to 
assess the chance that patches of high risk fuel arrangements might 
result in unacceptable fire behavior. 
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INTRODUCTION 


Fire behavior in a woody fuel array such as found in the forest is a complex 
physicochemical process that does not lend itself to a simple solution. The process 
is further complicated by the spatial nonuniformity of the fuels involved in the 
combustion process. A fire responds to fuel nonuniformity by changing its rate of 
spread and intensity. As a result there is a distribution in the rate of spread 
and intensity experienced by the fire as it spreads through the fuel array. 


For the purposes of this study nonuniform fire behavior is predicted by modeling 
fire spread through a hexagonal network of fuel cells. Fire spread is assumed to be 
a process of contagious growth between cells. Fuel properties are allowed to vary 
from cell to cell in a prescribed manner but have uniform properties within the 
cell. Consequently, the nonuniformity of the actual fuel array is simulated through 
cell to cell variations and has a resolution limited by the cell size. Because of 
the nature of the modeling process it is necessary to devise a scheme for collect- 
ing data describing nonuniformity and a scheme for filling the hexagonal cell array 
in a manner that simulates the actual fuel nonuniformity. 


Historically, it was necessary that uniformity be addressed first. This ap- 
proach laid the groundwork for present development. Rothermel (1972) chose a path 
of research into fire behavior that allowed the result to be applied to the needs 
of the fire manager. Rothermel's variables were the fuels and the environment in 
Which they were found. That same view exists here. The fundamental approach of 
the model described here is the assumption that small portions of a nonuniform fuel 
bed can be considered uniform. Nonuniform fire behavior then can be examined by 
following the progress of the fire as it changes speed while moving through differ- 
ent but uniform subunits of the area. This model was developed as part of an inte- 
grated effort for solving fire behavior problems at the Northern Forest Fire Labora- 
tory. The model is not intended for direct simulation of actual fire situations 
but rather is offered as a means to develop simplified methods for solving fire 
behavior problems in the field. 


Some aspects of fuel nonuniformity such as occasional absences of fuel have a 
direct effect on fire behavior. However, most often fuel nonuniformities are more 
subtle and must be viewed through another interpretive system, a model of fire be- 
havior. Consequently, we define fuel nonuniformity in terms of fire nonuniformity 
through a model of fire behavior that responds to spatial variations of the fuel 
array. 


Uniform fuel arrays and uniform fires do not occur. All fire behavior in 
woody fuel arrays is nonuniform. Fire behavior becomes uniform by definition. As 
a point of comparison, most people consider fire spreading over a sheet of paper 
as uniform. But, it may be viewed by some as nonuniform by focusing on the elemen- 
tal combustion process. Certainly, fuel uniformity does not exist in natural woody 
fuel arrays. The arrays are comprised of separated particles. Fons (1946) described 
fire spread as a series of ignitions. Fire spread appears to follow this description 
on the scale of the fuel separation. The flame ignites the particle and then spreads 
along that particle until it reaches another particle. The flame may ignite the next 
particle by flame contact or the particle may ignite after radiation from the flame 
front has preheated the fuel sufficiently to drive off combustible gases that burn 
when in contact with flame. 


A fire appears uniform when the particle separation is small compared to the 
flame size (or the range of influence of the flame). Individual flames coalesce 
into a solid moving front and the fire spreads in response to the bulk fuel proper- 
ties of the array--load, particle size distribution, and depth. These properties 
are defined within a volume element, averaging out the small variations of particle 
separation. Averaging is not limiting so long as the volume element encountered by 
the flame is small compared to the size or range of influence of the flame. It is 
this level of nonuniformity that we address. Wind and slope can act to increase the 
size of the flame and orient it so that the spatial variation of bulk fuel proper- 
ties is small compared to the area affected by the flame. Thus, a nonuniform fire 
can become uniform in the presence of wind or slope. 


Although resolution is limited by the cell size, the introduction of wind 
alters the overall perspective of resolution. Wind increases the size of the area 
sustaining active combustion and therefore increases the error of locating the 
fire. Consequently, in the presence of wind, the size of the cell can be increased 
without altering accuracy in locating the fire. 


The model is applied to two examples of nonuniformity: slash, residue left 
after tree harvesting, and a mixed community of grass and sagebrush. (Slash is 
more uniform than a mixed community of grass and brush.) The reader should expect 
to gain an appreciation for the problem of defining fire nonuniformity, developing 
a model of fire behavior that responds to fuel nonuniformity, and an appreciation 
of the kind and form of the results obtained from the nonuniformity model. Attempts 
are not made to validate the model. Consequently, the rigors of replication are 
replaced with a logical flow of fire behavior concepts--concepts derived from the 
uniform fire behavior model. The initial advantage is a consistent manner of 
handling field data that is applied to nonuniform fire behavior. 


The response of a spreading fire to the bulk properties of woody fuel arrays-- 
as found in forest fuels--has been investigated by Fons (1946), Thomas and Simms 
(1963), Anderson (1968, 1969), Frandsen (1971), Steward (1971), and Rothermel (1972). 
Rothermel incorporated fuel parameters (load, size, depth) and the fuel interactions 
into a model of fire spread through a continuous fuel array. Although the fuel 
array is continuous, it may be heterogeneous in size and type. Live and dead fuel 
may be included if mixed in the same stratum. 


Fuel parameters for the uniform fire spread model are categorized as living or 
dead, and averaged within a specific set of fuel size classes. The fuel array is 
assumed to be continuous. As a consequence, the model given by Rothermel (1972) 
predicts well for spatially continuous fuels and becomes increasingly less accurate 
as fuel discontinuity increases. To properly assess fire behavior, continuity must 
be included as an essential parameter in the mechanism of fire spread. Brown (1966) 
described the problem of continuity as follows: 


Closely related to compactness but on a larger scale is con- 
tinuity or patchiness of fuels. It represents the degree of 
change, horizontally and vertically, in the physical character- 
istics of fuels existing over a given area and is a measure of 
the uniformity of continuous combustion (for a constant set of 
weather conditions). At present, a meaningful quantitative de- 
scription of continuity is lacking. Development of an objective 
description of continuity would be of benefit to continued fire 
research as well as fuel appraisal and fire control planning. 


A discontinuous fuel array exhibits abrupt changes and is a unique example of 
nonuniformity. We employ a broad definition of nonuniformity that includes discon- 
tinuities. 


A simple example of a nonuniform array is one in which only the depth varies. 
The fire cannot achieve a constant rate of spread throughout the entire fuel complex, 
but may achieve it for some uniform subunit of the larger nonuniform array. Conse- 
quently, fire will accelerate and decelerate as it moves through the array. An ob- 
server usually calculates the average rate of spread from the time it takes the fire 
front to travel some given distance, but the result may not be related to either the 
individual depths or to the average depth. 


Nonuniform fire behavior implies more than one result for the rate of spread 
and intensity. Results should be expressed as a frequency distribution allowing the 
user more complete information on which to base a decision. The breadth of the 
distribution indicates the range of options to be considered in the management or 
control of a fire. Rothermel (1974) gives general accuracy requirements for the 
application of fire behavior models to fire management and control ranging from 
training aids to real-time fire predictions. As a training aid the new information 
will help emphasize the probabilistic nature of fire behavior. The impact on planning 
and management can be profound allowing a realistic assessment of the range of effects 
for alternative treatments of the land. The highest requirement for accuracy is pre- 
dicting real-time fire behavior. A knowledge of fire nonuniformity at this stage is 
essential. 


A method of collecting information from the field that represents fuel nonuni- 
formity is not common to present. fuel inventory systems. Use of average fuel para- 
meters in the uniform model as an alternative produces less reliable results as the 
fuel array becomes more nonuniform. An averaging of the fuel parameters prior to 
processing by a fire behavior model ignores the variable nature of fire as it moves 
through a nonuniform fuel array. A change in the fuel does not imply a proportional 
change in the fire behavior. An improved estimate of fire behavior can be derived 
from an analysis of the distribution of fire characteristics produced by the model 
as the fire passes through the array. 


Both rate of spread and intensity are implied in the frequency distributions of 
fire behavior results. Distributions of the rate of spread allow for a realistic 
assessment of the actual range of spread rates and area growth rates essential to 
prescribed burning and the control of wildfire. Distributions of the fire intensity 
provide insights for an assessment of the distributions of flame lengths and crown 
scorch heights to be expected on a site. As research into quantification of the 
heat pulse impinging on the site progresses and is related to the intensity distri- 
bution, assessment of the proportion of the burned area that will regenerate vege- 
tation after a fire should be possible. Variation in regeneration is seen in the 
response of seedlings in the vicinity of pile burning. Davis (1959) writes: 


It is common to observe that seedlings and other vegetation often 
do not become established for several years where large slash piles 
have been burned. A profusion of vegetation may fringe the area, but 
the center may be bare. 


These efforts also should provide inputs to hydrologic and soil stability impacts. 


SIMULATING FIRE BEHAVIOR 


We are concerned only with nonuniform fire spread in the horizontal plane. 
Arrays with large vertical nonuniformity are not considered. The results of the 
Simulation must reflect the variability in fire behavior due to the spatial non- 
uniformity of the fuel. Fire spread results are presented as a frequency distri- 
bution of the rate of spread of the fire front at variable lapsed times from igni- 
tion of the initial fire front. Intensity distributions are presented for those 
lapsed times. 


The simulated fuel bed consists of an array of cells. Each cell is described 
by a set of basic fuel parameters (Rothermel 1972) that affect the rate of fire 
spread. The time to move from cell to cell depends on the parameters of the two 
adjacent cells. The fuel is assumed to be evenly distributed within each cell. 


The fire begins as a line source and travels from cell to cell by contagious 
growth through a series of ignitions and spreads at a rate based on a minimum delay 
time (appendix I) within each individual cell (fig. 1). Hexagonal cells are used 
because they do not have point contacts when arranged in an array and offer the 
maximum number (six) of growth directions while maintaining a constant distance 
between cell centers. 


The delay time is the core of this analysis. Delay time is viewed as the time 
it takes for the fire to spread through the cell. Thus it is the amount of time the 
fire is delayed in a cell before it can attack an adjacent cell. The delay time is 
the sum of two parts: the time to achieve the quasi-steady spread rate in the cell 
and the time to spread the remaining distance at a steady rate of spread. The time 
to reach the quasi-steady state is assumed to be the residence time of the cell and 
the remaining time is the remaining distance divided by the quasi-steady rate of 
spread for that cell. Consequently, the cell chosen must be of a size such that the 
time for the fire to spread through the cell is equal to or greater than the resi- 
dence time of the fire in the cell. An essential restriction of the model is that 
the influence from a burning cell extend no further than its immediate neighbors. 


After passage of the front, the fire will remain burning in the cell for a 
period of time dependent on cell fuel properties, the properties of the cell it 
was ignited by, and the size of the fuel cell. The delay time can be shortened 
after ignition if heat from another adjacent cell has sufficient intensity to 
offer a delay time less than the present waiting time. 


Fire from one cell can ignite an adjacent cell after it has reached its es- 
cape level, i.e. waited a time equivalent to its delay time. If the adjacent cell 
is presently unignited, the fire moves in, delay time is assigned, and the fire be- 
gins its waiting period to reach its escape level. If the adjacent cell is already 
ignited and the proposed delay time is less than the present waiting time, then the 
waiting time is replaced by the delay time. Otherwise the cell retains its waiting 
time. A detailed discussion of the delay time appears in appendix l. 
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Figure 1.--Fire front progress through two nonuniform hexagonal fuel arrays. 


The rate of fire spread also has a dependence on wind and slope. Rothermel 
(1972) used only the maximum effect, i.e. wind and slope in the same direction. 
Provision was not made for combining nonparallel wind and slope nor for the compu- 
tation of rates of spread in the six spread directions required by the hexagonal 
array. Albini, ! however, developed a method for combining wind and slope that 
gave the resultant magnitude of rate of spread as a function of an angle relative 
to the directions of the wind and slope. With Albini's model, it was possible to 
compute potential rates of spread in all six directions in the presence of wind 
and slope. 


1Ajbini, Frank A. Memorandum, subject, Combining wind and slope effects on 
spread rate, to R. C. Rothermel, Northern Forest Fire Laboratory, Missoula, 
Montana 59806, Jan. 19, 1976. 


As a comparison, the Dijkstra (1959) algorithm employed in the Kourtz-0O'Regan 
(1971) fire spread model locates the paths of least resistance based on delay times 
dependent only on the fuel descriptors of each individual cell, without dependence 
on an adjacent cell. As each cell is reached--ignited--by the fire, the lapsed time 
Since fire starting time is recorded in the cell. Isochrones can then be drawn to 
illustrate the progress of the fire. The two models are similar except for the 
differences in delay time mentioned above. The Dijkstra model does not require up- 
dating since it maps the path of least resistance to fire spread. Presently the 
Kourtz-O'Regan model uses very large cells and average fuel parameters to describe 
the average rate of spread and thus the delay time to consume the fuel cell. The 
hexagonal model operating on cells that are subunits of the larger Kourtz-O'Regan 
cell then can provide a distribution of rate of spread values needed to compute 
the probable time a fire takes to consume a cell in the Kourtz-O'Regan model. Thus, 
data obtained from the hexagonal model could be used as input to the Kourtz-O'Regan 
model. 


For a continuous fuel array the fireline intensity is the product of the reac- 
tion intensity and the combustion zone depth (Albini 1976) and assumes a constant 
reaction intensity throughout the combustion zone. For the nonuniform array the 
reaction intensity is assumed constant throughout the cell but may vary from cell 
toy.ceulie 


Following a suggestion by Frank Albini of the Northern Forest Fire Laboratory, 
the fireline intensity is obtained by summing the products of the reaction intensity 
and that portion of the cell contributing to the combustion process. Each column is 
scanned perpendicular to the initial fire front (fig. 1): 


Le » Co) ae 


where I, is the column fireline intensity, (1,); 
ie cell, FS (appendix II) is that fraction of the es cell that is contributing to 


is the reaction intensity of the 


the combustion process, and D is the cell width. All column intensities then are 
grouped to form a frequency distribution. 


The sum of the products, SS F.D, is the combustion zone depth if the fire is 
AE ’ 
burning perpendicular to the initial line of fire. Occasionally, portions of the 
fire front may be burning to the side. Scanning down the columns would then give 
some combustion zone depths that are unreasonably high. This results in spikes 
in the distribution at high fireline intensities that can be easily located and 
disregarded. 


Fuel Array Assessment 


The geometric nature of the fire behavior model presented resulted in a fire 
spread algorithm coupled to a hexagonal array of fuel cells. However, application 
to the field is not possible unless we are able to fill the hexagonal array in a 
manner that preserves the horizontal stratification of the actual fuel array. 
Filling the hexagonal cells requires fuel array data (load, size, depth) that must 
be acquired in a manner that fills the needs of the cell filling algorithm. In 
general, we should look for some classifiable character to give the assessed non- 
uniformity a recognizable distinction related to its habitat type or other compar- 
able classification, i.e., fuel type, and age. Habitat classification according 
to Daubenmire and Daubenmire (1968) is presently used by the USDA Forest Service 
in Idaho and Montana to classify vegetation and its associated environment. Fuels 
are a byproduct of the habitat type but may occur in different arrangements of 
load, particle size distribution, and depth, and thus are classified separately. 
Fuel models used in the U.S. National Fire Danger Rating System are obtained by 
grouping depth, particle size, and load into a classification scheme. Data 
gathered from sampling fuel arrays should characterize the horizontal patter 
of differing fuel types as well as the spatial occurrence of the basic fuel pro- 
perties (load, size, depth) within a fuel type. Other fuel properties--heat con- 
tent, mineral content, and particle density--have little variation within the fuel 
type and are assumed constant. Fuel moisture is time-dependent responding to diur- 
nal changes in humidity and temperature. However, these changes are sufficiently 
slow so that fuel moisture can be considered constant over the duration of the 
fire being examined. Although moisture can be introduced as a spatial variable, 
it is held constant so that the response to nonmoisture fuel variability can be 
emphasized. It is anticipated that experienced land managers will make adjustments 
for moisture changes. 


Two methods of fuel assessment are presented; (1) sampling a specific size 
area at periodic intervals along transects, and (2) if components are random, 
evaluating the percent cover and describing the uniform fuel properties of each 
component. 


SLASH 


The first method stated was employed in slash owing to the absence of recog- 
nizable patterns in the spatial arrangement of the fuels. Transects were obtained 
from a slash area composed primarily of western larch and grand fir. The area was 
essentially a clearcut with only a few remaining trees. Trees were cut to an 
8-inch (20 cm) unmerchantable top and the entire tree except the top skidded to 
the landing. Nonuniformity was assessed in terms of load and depth along a 100-foot 
(30.5 m) transect at 2-foot (0.61 m) intervals. The fuel load was estimated by size 
class from the number of intercepts through a vertical sampling plane (Brown 1974). 
The following size classes were assessed: lh, 10h, and 100h.2 Pieces greater than 
3 inches in diameter were measured but not considered in the model because they do 
not significantly contribute to fire spread. 


*Fuel size classes are characterized by the time lag constant related to their 
ability to respond to humidity by absorbing or desorbing moisture (Fosberg 1970). 
The 0-1/4 inch (0-0.63 cm) size class is called lh, 1/4-1 inch (0.64-2.54 cm) the 
10h, and the 1-3 inch (2.55-7.62 cm) the 100h. 


Because of the roughness of the surface of slash, it was necessary to define 
a depth called the bulk depth. It is an estimated mean that retains the bulk den- 
sity of the inventoried fuel load and is confined to a cylinder 1 foot (0.30 m) in 
radius whose central axis is perpendicular to the slope at the sample point. Four 
estimates are made within the cylinder, each representative of 25 percent of the 
area (not necessarily a quadrant). Vertical gaps of more than 1 foot are subtracted 
before each estimate is completed. Gaps less than 1 foot are assumed to maintain 
vertical continuity through the potential flame height. The average of the four 
estimates is recorded as the estimated mean depth. 


Each logged unit was measured in two stages. A grid of approximately 30 sample 
plots was established for each unit. In the first stage, the load for each of the 
four size classes was measured at each plot according to Brown (1974). In the 
second stage, the unit was arbitrarily divided into four sections and a sample plot 
was chosen at random from each section. At the first sample plot, a 100-foot 
(30.5 m) transect was selected at random and oriented in one of three directions: 

(1) 60° counterclockwise to the uphill direction, (2) uphill, and (3) 60° clockwise 
from the uphill direction. A random selection of the two remaining orientations 

was made at the second plot. The last remaining orientation was made at the third 
plot. Orientation on the fourth plot was chosen at random from all three directions 
as in the first plot. Sample planes 2 feet (0.6 m) across were oriented first along 
the transect line and then perpendicular to this line with the sample point as the 
center. The bulk depth was recorded for each sample point (fig. 2). The 10h size 
class load was measured at both orientations of the sample plane at one-third of 

the sample points which were selected at random. The range of depths and loads ob- 
tained from transects are given in figures 3 and 4. 


In stage two, only the 10h intercepts were counted. Data gathered in the first 
Stage provide the information for relating the load in the 10h class to the other 
two classes (lh and 100h). 


The resulting 2-foot sampling interval established a minimum cell size for the 
hexagonal array. A smaller size would seriously degrade the accuracy of the load 
estimate obtained from the planar intersect technique. If larger cells are re- 
quired, they can be constructed by combining cells. 


Three of the four transects taken on the slash area had similar cumulative 
depth distributions. The remaining transect was discarded leaving a data base 
of 150 bulk depths and approximately 50 10h-fuel load estimates. The unit was 
generally described as having a light slash load. Consequently, this unit 
should exhibit correspondingly low rates of spread and intensities. 
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Figure 2.--A representative transect of the bulk depth evaluated every 2 feet. Three 
orientations are shown: (1) 60° coumterclockwise to the uphtll direction, (2) up- 
hill, and (8) 60° clockwise to the upntll direction. 
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Figure 3.--Frequency distribution of the slash bulk depth. (All measurements were 
made at a resolution of 1 tneh (2.54 em). 
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GRASS AND SAGEBRUSH 


The second method, evaluating percent cover, was used for the mixture of grass 
and sagebrush because this mixture can be separated according to load breakdown by 
size class and depth, resulting in two components. Grass is predominantly made up 
of single sized particles. The bulk density, load divided by depth, varies from 
0.028 to 0.085 lb/ft? (0.45 to 1.36 kg/m?) with a mean of 0.05 lb/ft? (0.80 kg/m?) 
according to data Sneeuwjagt (1974) collected from the Soil Conservation Service 
on Western United States grasses. A knowledge of either the load or the depth is 
sufficient to quantify the amount and arrangement of the grass fuel if the bulk 
density is known. For this example, we have considered all of the grass fuel to 
be dead. 


Studies by Rittenhouse and Sneva (1977) and Brown (1976) allowed the evalua- 
tion of fuel loads for sagebrush according to height and largest planform diameter. 
The first study related the mass to the crown dimensions of the plant. The second 
study related the mass of single stems (including branches and foliage) to the 
basal diameter of the stem. But the second study (Brown 1976) went a step further 
to provide a breakdown of the mass by size class. Assuming that the sagebrush 
plant originated from a single stem. We were able to compare the two methods and 
OD een Ae: load (mass per planform area) related to the crown dimensions of the 
plant. 


The following dimensions were chosen to represent sagebrush; an average height 
of 3 feet (91 cm) and an average diameter of 2.5 feet (76 cm). Table 1 lists the 
fuel parameters associated with these bulk dimensions. Fuels larger than 1/4 inch 
(0.64 cm) have been ignored in this analysis because they have only a minor impact 
on the spread rate and fireline intensity; these larger fuels are not consumed in 
the initial combustion process and therefore do not contribute to the spread process. 


Table 1.--Fuel parameters for grass and sagebrush 


: Surface/volume : Dry load * Moisture! * Moisture of2 


Fuel description Be 3 gine IB Ees : k ina : content extinction 


Grass (fuel depth = 1.00 ft (0.30 m)) 
3000 Vex OO275 Ons: 0.05 ORS 


Sagebrush (fuel depth = 1.48 ft (0.45 m)) 


Li re foliage 1500 49.2 0.0538 0.263 1.00 2.00 
Lave lh 677 LONE: 0.0941 0.459 0.50 2.00 
Dead lh 677 DD ie. 0.0235 0.115 0.09 0.20 


The values of fuel parameters common to all fuels discussed are: 


Low heat vadiennns imam.) mer bc as | eae SOOOKORBRUAI DIN Gls 595) Jeii/icg)) 
Particle idensaty, au.) Sct Ot neers uted ar 32.0 lb/ft? (0.51 g/cm?) 
Fractional | mineral angen Se ais Spon ote 0.06 

Fractional! effective mineral content: bi tolls ai ere 0.01 


lFraction of dryweight. 
2A parameter involved in the computation of the moisture damping coefficient 
(Rothermel 1972). Higher values allow the fire to spread at higher moisture contents 


3Data on file at the Northern Forest Fire Laboratory, Missoula, Montana. 
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Data collected and analyzed? according to Holgate (1965) shows sagebrush plants 
to have a random distribution. Consequently, the additional information required to 
fill the hexagonal fuel array is the number density of plants. The density, coupled 


with the average planform area, dictates the percent cover that must be achieved for 
a random filling of the array. 


Fuel Array Generation 


SLASH 


Assignment of the bulk depth to each 2-foot hexagonal cell is the first step 
toward generating the slash fuel array. The fuel load then can be assigned to each 
cell through direct and indirect relationships to the depth (see appendix III). 


After surrounding the area to be simulated by a boundary of depth values 
having the same distribution as the array but not the spatial order, the remain- 
der of the area can be filled with depths using a linear model suggested by Mike 
Marsden of the Northern Forest Fire Laboratory: 


en De =P bX + b,x, + bX. +e 


where Y is the depth being assigned to a cell and X,, X,, and X, are the depths that 
have been assigned to three adjacent cells and e is the error. 


(n 
CORD) 
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The dependent variable, Y, and independent variables, X's, change as the array is 
filled. A more extensive array using the same principle is shown below. 


M4 X, X, X, 
Yy a b c 
Y5 c d e 
Ys e £ g 
Y4 g h a 
Ys 71 5 2 
Y6 Y5 € Y3 


The coefficients bo bi b,, and be are derived from the mean bulk depth and 
Serial correlations of lags 1 and 2 of the bulk depths along linear transects 
through the slash array (appendix IV) and not according to the cell filling for- 
mat, shown above. The following data were obtained from field transects and are 
used to develop the array for the slash example. 


Average depth = 4.44 in (11.3 cm) 


Variance = 23.01 ine (148.5 cm?) 


Serial correlations 


Orientation Lag 1 Lag 2 
i 0.45 -0.044 
2 0.69 -- 
5 0.62 “Oo eke) 


Serial correlations are simply correlations of data pairs (Snedecor and 
Cochran 1967) obtained from sequential transect depths of lag 1 and 2. A lag 
of 1 designates a correlation of adjacent depths, whereas a lag of 2 designates 
a correlation of data pairs obtained by skipping one depth in the transect se- 
quence. The numerical value of the correlation is the correlation coefficient 
that allows the user to quantify the similarity or dissimilarity of these data 
pairs. 


The error, e, associated with each prediction of y is obtained through ran- 
dom access of the cumulative distribution of the bulk depth (see appendix V). 
Because of the error term, generated fuel arrays are not identical: This 1s in 
agreement with the goal to produce a pattern that has the essence of the array 
but is not an exact reproduction. An example is wallpaper design. Your eye 
recognizes a pattern, but may not find exact comparison. 
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The average fuel parameters of the simulated fuel array are given in table 2. 
A comparison with figures 3 and 4 shows the tabular load and depth values to be 
within one standard deviation of the averages given in the figures except for 100h 
fuels. The foliage load is not given as a distribution. The foliage load is a 
fraction of the sum of the lh, 10h, and 100h fuel loads (appendix III). The simu- 
lated serial correlations of the depth for lag 1 were 0.12, 0.31, and 0.11 for 
orientations 1, 2, and 3, respectively. These correlations are reduced consider- 
ably from the original correlations used to generate the fuel array. However, 
these correlation data showed the best comparison while maintaining an average 
depth comparable to the original data. Other correlations can be obtained by 
manipulation of the error function, but this in turn causes unacceptable changes 
in the average depth and its variance. The average depth of the final simulated 
array is 5.45 inches (13.8 cm) with a variance of 18.03 in* (116.3 cm*). These 
data are in good agreement with the original depth data used to generate the fuel 
array. 


Table 2.--Average fuel parameters for slash 


— . Fe «- 
: Surface/volume ; Dry load ~: Moisture! 


Fuel description: yee el Sa 1b/£t? : ken : content 


Slash (fuel depth = 0.43 ft (0.13 m)) 


Dead foliage 2000 65.6 0.118 0.576 0.05 
Dead lh 436 14.3 0.0363 Ooh 0.05 
Dead 10h OI 3.0 0.246 1.199 0.05 
Dead 100h 29 LAG 0.190 0.926 0.05 


Moisture of extinction = 0.25 


lFraction of dryweight. 


GRASS AND SAGEBRUSH 


Generation of the grass and sagebrush fuel array is less complicated than 
generation of the slash fuel array because of the random distribution of the sage- 
brush plants within the grass matrix. The cell size of the fuel array is chosen 
to approximate the average diameter of the sagebrush plant, 2.5 feet (76 cm). It 
is necessary then only to classify the cells randomly as grass or sagebrush from 
a distribution that reflects the percent cover of the sagebrush within the grass 
matrix. A value of 30 percent was chosen as representative of the percent cover 
of sagebrush. Cells were assigned a fuel type of either grass or sagebrush. The 
fuel descriptors are given in table 1. The fuel depth listed for sagebrush is an 
equivalent height that when divided into the fuel load gives a bulk density equiva- 
lent to the shrub crown. 
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RESULTS 


Predictions of the fire spread rate and intensity were obtained for slash at 
0 and 2 mi/h (0 and 3.2 km/h) and for the grass-sagebrush mixture in the absence 
of wind. Slope was 0 for both examples and the wind was perpendicular to the ini- 
tial line source. The important difference from other forms of the result is that 
the predictions are presented as distributions. 


The rate of spread distribution was obtained from the distance traveled in 
a specified time along each column of cells perpendicular to the initial fire 
front (fig. 1). Furthermore, the characteristic distance that the fire has 
traveled through a cell--and thus the fraction consumed--can be calculated at a 
given time (appendix II). Consequently, the accuracy of the overall distance 
traveled is not limited to the cell size. 


The predicted intensity, Byram's fireline intensity, is presented as a distri- 
bution made up of the intensities from each column of the array. 


SLASH 


The distributions of the rates of spread at 0 windspeed and at a 2 mi/h wind- 
speed and 0 slope are given in figure 5. A prediction of the rate of spread assuming 
a uniform fuel array having the average fuel parameters of the simulated array as 
given in table 2 is presented also for comparison to the distribution. The average 
spread rate of the uniform model lies 31 percent” below the nonuniform model (average 
of the distribution) at a 0 mi/h windspeed and 21 percent below at a windspeed of 
2 mi/h. The overall range for both cases lies from 65 percent below to 76 percent 
above the nonuniform model average. 


“The nonuniform model average was used as the base for all percentage calcula- 
tions. 


U5 


Wind, 2 mi/h (3.2 km/h) 


20 


5 
> 
oO 
= 
am UNIFORM MODEL 
g 10 (FUELS AVERAGED) 
i NONUNIFORM 
MODEL AVERAGE 
5 
0 1.0 2.0 3.0 
RATE OF SPREAD (ft/min) 
(m/min) 
0 0.1 0,2 0.3 
20 
No Wind 
15 
> UNIFORM MODEL 
2 (FUELS AVERAGED) 
S 10 
re NONUNIFORM 
= MODEL AVERAGE 
5 
0 0.2 0.4 0.6 0.8 1.0 1.2 


RATE OF SPREAD (ft/min) 


Figure §.--Frequency distribution of the rate of spread in slash fuels at 0 wind- 
speed and at a 2 mt/h_windspeed. The average rate of spread is tndicated for 
the wut form model (average fuel parameters from the stmulated array) for com- 
partson with the nonuniform model (average value of the distribution). 
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The distributions of the fireline intensities at 0 windspeed and at a 2 mi/h 
windspeed and 0 slope are given in figure 6. A prediction of the intensity from 
uniform fuel properties is again presented for comparison, similar to the presenta- 
tion of the rate of spread prediction. The average intensity of the uniform model 
lies 43 percent above the nonuniform model at 0 mi/h windspeed and 29 percent above 
at 2 mi/h. Both distributions decrease monotonically with increasing intensity and 
are bounded on the low side by 0 intensity. The upper limit is 6 times the nonuni- 
formity model average at 0 mi/h windspeed and 7 times that average at 2 mi/h. 
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FIRELINE INTENSITY (Btu/s/ft) 


Figure 6.--Frequency distribution of fireline intensity in slash fuels (0 windspeed 
and 2 mt/h windspeed). 


Wy 


GRASS AND SAGEBRUSH 


The distribution of the rate of spread in the absence of wind and slope is 
given in figure 7. Uniform model predictions are obtained by averaging the fuels 
according to the percent area covered in the fuel array. In addition to the uni- 
form model average and the nonuniform model average, the uniform rates of spread 
for grass and sagebrush are given individually along with their average weighted 
by percent area covered. The additional predictions are possible because the 
fuel array is made up of two distinct components. 


In addition to the uniform and nonuniform model averages as presented 
in figures 5 and 6, the uniform results for grass and sagebrush are 
indicated separately along with the average, weighted by percent area 
covered. 
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Figure 7.--Frequency distributions of rate of spread and fireline intensity in a 70 
percent grass-30 percent sagebrush mixture (tn absence of wind). 
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The rate of spread distribution lies in a narrow band near 3 ft/min. See 
figure 7. All values presented are less than 4 ft/min. All three averages-- 
uniform, nonuniform, and weighted--lie within the extremes of uniform sagebrush 
and uniform grass. 


The fireline intensity distribution extends from 0 to 30 BTU/s/ft and in- 
cludes the uniform model averages. In contrast to the rate of spread, both the 
uniform and nonuniform model averages lie outside the extremes of uniform sage- 
brush and uniform grass--exceeding both values. 


Fires were not simulated in the grass-sagebrush mixture at 2 mi/h windspeeds 
because the transition zone depth exceeded the cell size, 2.5 ft (76 cm). See 
appendix I. The grouping of cells into larger cells is possible but in this case 
the fundamental cell is the size of a sagebrush plant. Averaging such dissimilar 
fuel cells at this time does not serve the purpose of this paper--to illustrate 
the nonuniform modeling approach and the form of the result. 


DISCUSSION 


Fuel nonuniformities, and resulting nonuniform fire behavior, although recog- 
nized as a problem, have been ignored because of a general inability to handle non- 
uniformity both in its assessment and simulation and in the form of the result. His- 
torically we have attributed a single rate of spread to a unit of land. This concept 
existed before the formulation of the uniform fire spread model (Rothermel 1972). 

The assumption of uniformity in the first attempts to describe firespread in wild- 
land fuels was a logical first step. The assumption may have enforced our view of 
Single valued results but is not the source of that viewpoint. The recognition of 
fuel nonuniformity carries with it the realization that fire exhibits nonuniform 
behavior. Consequently the result must be in the form of distributions, excursions 
about the average. 


Two fuel arrays, slash and a grass-sagebrush mixture, are presented as examples 
for describing the nonuniform fire behavior model. Slash is a single fuel type with 
varying depth and load. The grass-sagebrush mixture is a mixture of fuel types. 
Unlike slash, most any location in the mixture has a dichotomous description as 
either grass or sagebrush. 


All possible combinations of obtaining average fire behavior results using the 
uniform fire model were employed for comparison with the distributions. For slash, 
the average depth and fuel load were used. For the grass-sagebrush mixture, there 
were the additional results obtained from uniform grass, uniform sagebrush, 
weighted average of grass and sagebrush fuels, and the weighted average of the 
two results for uniform grass and uniform sagebrush. 
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In general, there is agreement among the averages of the uniform model (fuel 
parameters averaged) and the nonuniform model (distribution averaged). The trend 
indicates that the uniform model average is less than the nonuniform model for 
rate of spread but greater for fireline intensity. The higher frequency of zero 
fireline intensity in slash is attributed to the frequency of cells not having 
fuel. This is consistent with the grass-sagebrush mixture results where empty 
cells were not present. 


The average fireline intensities of the slash are low values that suggest a 
minimum capability to spread fire. However, the distribution indicates that por- 
tions of the fire at a 2 mi/h windspeed would be comparable to a prescribed litter 
fire but would be easy to control. The last statement illustrates the importance 
of the distribution. Although the results are not alarming, a single average value 
even if correct does not indicate the upper limits of the intensity. 


The upper limit on windspeed can be overcome by averaging cells into larger 
cells so that the influence of the fire from one cell does not extend beyond adja- 
cent cells. This is in keeping with the relative reduction in the influence of 
Spatial variations in fuel as the flame increases in response to wind. Future 
applications will likely employ general fuel models having flexibility in resolu- 
tion as needs arise.° 


The form of the result offers new ways of presenting alternatives to the land 
manager. With a distribution of results it is possible to ask what the probability 
is that a given range of rate of spread or intensity will occur. A high rate of 
Spread or intensity although having a low probability may be intolerable. Viewing 
fire behavior in the form of probabilities gives a more accurate description of 
the probable fire behavior and a reasonable basis for management decisions. 
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NOMENCLATURE 


Definition 


Reaction intensity 
Fireline intensity 
Burning fraction of cell 
Cell width 


Transition distance where fire is influenced 


by old and new cell 


Rate of spread; or, 


Multiple correlation coefficient in statistical 


analysis location 
Delay time 
Residence time 


Surface area-to-volume ratio 


o value within cell weighted by total exposed 


surface area of each particle size 
Particle diameter 
Time 
Rate of spread in transition distance 
Effective heating number 
Volumetric heat of preignition 
Propagating heat flux ratio 
Wind factor 
Slope factor 
Distance fire has traveled in cell 


Time dependent decreasing influence from 
old cell 


Time dependent increasing influence from 
new cell 


Time since ignition of cell 


Error term 
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none 


min 


none 


APPENDIX I: DELAY TIME 


The delay time of a cell is the time that must elapse from the time that the 
cell is ignited until it is able to ignite adjacent cells. In the simulation, we 
are concerned with fire behavior only at the resolution of the cell size. The 
fire is viewed as jumping into a cell at ignition and then waiting out the delay 
time before jumping into unignited adjacent cells. However, to develop a method 
of calculating delay time, the behavior of the fire within a cell must be considered. 


After cell (Gi+l) is 1gnited by cell 1, the fire spreads a transition distance, 
D', before it reaches a quasi-steady state. During the time that it takes the fire 
to spread this distance, the influence of cell i is decreasing and the influence of 
cell (itl) is increasing. The fire spreads the remainder of the distance through 
(i+1) at a rate which is calculated according to the uniform fire spread model 
(Rothermel 1972). 


The delay time is expressed as the sum of the residence time and the time 
that it takes the fire to spread the remaining distance (D-D') at a uniform rate: 


= = ! 
ta real AG D ORS gh 
where 
CARA = residence of cell (i+1) 
D = cell width 
D' = transition distance 
Rea = quasi-steady state rate of spread in cell (itl). 


The residence time, (ee 


point as it spreads through a fuel array. Residence time is assumed to be a mea- 
sure of the time that it takes a fire in cell (itl) to achieve the quasi-steady 
State where it is no longer influenced by cell i and spreads at a rate dependent 
only on the fuel parameters of cell (i+1). Our goal in constructing the hexagonal 
fuel array is to choose the cell size such that the time for the fire to spread 
through the cell is greater than the characteristic residence time of the fuel 
particles in cell (i+l). 


is the amount of time that fire exists at a given 


The residence time is expressed as the ratio of the combustion zone depth (ap- 
proximately the horizontal region of active flaming) to the rate of spread. Ander- 
son (1969) found the following approximation for the residence time in terms of 
the particle diameter d: 


eS Goll 
r 


where es is in minutes and d is in cm. The diameter, d, can be expressed in terms 


of the surface area-to-volume ratio: 


al a AY/Koj 
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Therefore, residence time can be calculated from the characteristic surface area- 
to-volume ratio: 


Cans 12.6/o 
NGS G es ne umes of Gm 
Each hexagonal cell has a characteristic surface area-to-volume ratio obtained 
from the average particle diameter weighted by the total exposed surface area of 
each size class of fuels within the cell (Rothermel 1972). Thus, the residence 
time is an average property of the cell. 


The values required for each cell for the delay time calculation are I> eli g? 
ae dy > and ‘a These values are calculated in turn from measurable fuel proper- 


ties. The reaction intensity, I, may be thought of as the "heat source" affecting 


Re 
adjacent cells. The second term, ea; o> is the product of the effective heating 


number, e€, and the volumetric heat of preignition, 4; . This term may be thought 


of as the "heat sink,'' the absorption of heat required for the fire to advance into 
the cell. The propagating flux ratio, €, is the heat coupling coefficient operating 
on the reaction intensity to obtain the propagating flux. The wind and slope fac- 
tors, o. and 5» are used along with the direction of spread in relation to the 


directions of the wind and slope to calculate a wind-slope factor, g(4 54.) -6 Cal- 


culation of I os 5 Ee d° and oe is described in Rothermel (1972). 
The rate of spread in cell (it+l) is initially influenced by the cell i. The 

fire is assumed to have achieved a quasi-steady state when the residence time of 

cell (i+1l) has elapsed. To smooth the change in rate of spread within cell (it+l), 

a gradual change is assumed. Two possible cases must be considered: 

i; Go Ga, adil (GS Gere hor icase Asthe einitlinencemotasceslitarr 

Tea see abl real re atsral 
terminates at time (tJ;- The influence of cell (i+l) continues increasing until 


time Cay) when the fire reaches its quasi-steady state. For case II the influ- 


i+l 


ence of cell i is terminated at time Cae rather than Caner At (Ce) the 


ard ataesl 

rear of the combustion zone leaves the boundary between the cells producing a 
burned out area between the cells whereupon cell i is assumed to no longer influ- 
ence cell (i+1l). Frank Albini of the Northern Forest Fire Laboratory suggested 
the following mathematical model for expressing these influences. 


The distance, s, that the fire has traveled in cell (i+l) by time t is: 


Case I: Gar SG) 


ie 1c 
5 = mJ f, (t)dt + Ria J meer (2) 


SAlbini, see footnote l. 
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Case II: (a > Cra 
t ic 
S = Ry f f,(t)dt + Steal h(t) dt (3) 
0 0 
where 
Lo t/ @_). f ORS te en hie) 
i = i 
f(t) = i (4) 
0 SEES Te: © (ar) 
ne aL 
iL es eH (Gee cope Te Se (He Do 
e) Shee hae (5) 
0 nope {ee Coa A 
Ce Dasa 10g (ES (han 
h(t) = (6) 
1 IEOHO 16S Geen 
t = time since ignition of cell (i+1) 
Ga, = residence time of cell i 
(Dara = reSidence time of cell (i+1) 


Re = rate of spread as the fire passes from cell i to cell (i+l) 
R. = quasi-steady state rate of spread in cell (i+l). 
The terms f(t) and f(t) are time dependent decreasing influences from cell i 
while h(t) is a time dependent increasing influence from cell (i+l). 


A graphical representation of the "influence" factors f, (t), f,(t), and h(t) 
are given in figure 8. 


1 1 
S f, (t) h (t) f,(t) h (t) 
— 
Oo 
te 
oO 
a 
= 
iL 
c= 

0 0 

(ae) (Fie) (Teale (T _) t 
rT eb S| TIME pee I ei d 


Figure 8.--Influencing factors governing fire transitton from cell i to cell (i+1) 


The quasi-steady rate of spread for cell (i+l) is calculated as described by 
Rothermel (1972): 


Rey = Cepia Sin G2 Ps) bier! OG iad: (7) 
The transition rate of spread as the fire moves from cell i into cell (i+l) is 
calculated from the heat source terms of cell i and the absorption terms of 
cell (i+1): 
This can also be written as: 
Rp Osea, (8) 


where 
C= (et eats, ea. eel 


Substituting equations (7) and (8) into equations (2) and (3) and evaluating 
the integrals over the intervals illustrated in figure 8, the distance that the 
fire spreads into a cell in a given time is as follows: 


Case I: (ree) ecea Gom 


) 


i+] 


CR(t i> 62/2 (wa) acean(E2/AGe oie e a fom 10) gees (9a) 
=) \CR(n),/2 Ree en eer ae Se Seen (9b) 
Die (t= Gere ik forse a: Fare aity (9c) 
Case II: Cee > Ce 
CR(t 22/2 Joo) 2 Rite Neato ere Cone (10a) 
S = 
Daca Coe aN Hor (cag < © SG, (10b) 
where R = R. 


ait 


The transition distance, D', can be obtained by solving equations (9b) and 


(10a) at t = (Cee 


Case I: (tJ, < (1 


Tee 
1 — alts: 


DS NEG © Qa? (11) 


26 


Case II: (t); > (Gy) 


15 4 
i a0 alent 


DS Ge) sg (Gah yc (12) 


Equations (9) and (10) can be solved for t to obtain the time that it takes 
the fire to spread a distance, s, in the cell. 


Case I: (ty SiGe Ea) ies 
Lc + Vco + 2B(1 - CEmN A C// Wkrae,. > Ole Daly (13a) 
for (ea) < (Oe OngGaxiel 
andy SJ RU) HE (t), / De, 2 
s/R, for Ca = (Casa (Geil (13b) 


and Q<s Ns) 5 


N25(CSL ATR - ss se Ce) A rare (13c) 


FOR ERMae)[C © (Cab) (ad ayaa //2 Gousiase! 


GED Cae (134) 
1059 MD SS S10) 
Case II: Ge > Cer 
2 
ae [-c + Ge sre A S(Ol 3] Au - OVE Ico Oe MO < Sox jy (14a) 
(s - D')/R + Cis LOT <Sa <i) (14b) 
where 
te i+] 
B= S/R(t 544° 


If the fire reaches its quasi-steady state in a cell, equation (13c) or (14b) 
for s = D is equivalent to finding the delay time as described in equation (1). 


As stated earlier, the cell size should be large enough for the fire to 
achieve a quasi-steady state in every cell. Otherwise, the influence of cell i 
would extend beyond cell (i+l). At present the algorithm does not handle this 
Situation. Increasing the cell size excessively may result in averaging out the 
nonuniformities we originally intended to examine. It is essential to design the 
array so that a high percentage of the hexagonal cells achieve a quasi-steady 
State rate of spread and accept the small error resulting from the few exceptions. 
If the calculated transition distance is larger than the cell size (D' > D), the 
delay time is calculated as the time that it takes the fire to travel the distance D. 


BI 


The following is a general procedure for obtaining the delay time. If 
Coa SE IGGE)) FADD substitute D for s, select the equation having the correct 


limits £rom equations 3a; b, om ©; and) solve forest. vite ( )2 = Gore. proceed 
: ‘ 5 re at ae atspll 
in the same manner using equation 14a or b. 


APPENDIX I: CELL BURNING FRACTION 


The fraction of a cell that is burning at time 6 since its ignition is evalu- 
ated from the distance equations for s in appendix I. Initially the front of the 
burning zone spreads into the cell as described by s with t = 6. The rear of the 
burning zone follows in the same manner but with t = 6 - res where = is the 


residence time in the cell. During passage of the fire, the front and the rear of 
the burning zone may exist in the cell together, separately, or be absent. The 
following equations account for all cases: 


F = 0 for 6<0 The cell has not been 
ignited. 
F = s(@)/D for D'<D and 0<6<t The fire front is in the 
a cell. There is no 
or for D'>D and 0<6dt , burned out area. 
F = (s (6) -s (6-1 ))/D for D'<D and t_<6<t Both the fire front and 
r — eid ; 
rear are in the cell. 
E = 4 for D'>D and ESC, Thement are celiieas 
burning. 
F = 1-s(6-t_)/D for D'<D and t ,<6<t ,+T The fire front has passed 
a6 = Gl sel" ae 
the cell. 
or for D'>D and SONG a Hur There is a burned Out area. 
F = 0 for OAS The cell is burned out. 
where 


F = fraction of the cell that is burning at time 8 
D = cell size 
D' = transition distance 


t, = delay time 


d 
ian cell residence time 
s(6) = distance the front of the burning zone has spread at time t = 6 according 


to equation 9 or 10 


(6-1, ) = distance the rear of the burning zone has spread at time t = oe 
according to equation 9 or 10. 
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APPENDIX Il: CELL LOAD EVALUATION 


The 10h load is assigned through a relationship between the load and depth. 
The lh and 100h are then assigned through two separate relationships of each of 
these loads to the 10h load derived from the first stage inventory. The necessary 
relationships are expressed in figure 9 as cumulative distributions of the follow- 
ing ratios: (upper) 10h load to bulk depth, (middle) lh load to 10h, and (lower) 
100h load to 10h load. For each cell in the depth array, the 10h load can be ob- 
tained by random access of the upper distribution in figure 9. The other two loads 
are determined in the same manner using the 10h load as a base and accessing the 
middle and lower distributions. 


An estimate of the foliage load was made based on the dominant species compo- 
Sition of slash, foliage retention by these species, and a knowledge of the foliage 
load relative to the sum of the lh, 10h, and 100h slash load. 


Dominant species composition of the slash was western larch and grand fir. Be- 
cause larch loses its needles quickly compared to grand fir and the slash had gone 
through one winter, it is reasonable to assume that western larch had lost all of 
its foliage while grand fir retained its foliage. Brown (1978) found that the 
grand fir foliage load was approximately 50 percent of the sum of the lh, 10h, and 
100h loads. Depending upon the relative amounts of western larch and grand fir, 
the foliage load could vary from zero to 50 percent of the overall sum of the lh, 
10h, and 100h slash load. An average of 25 percent was chosen. After lh, 10h, 
and 100h loads are determined for each cell, 25 percent of the sum is used to repre- 
sent the foliage load. 


APPENDIX IV: EVALUATING by, bi, bz AND bs 


Evaluation of the coefficients by Derry 10) 


and b. of the multiple linear 
regression equation: 


1 Be 


Yor bo + bX) + b,x, + b AX. +e 


is obtained by using serial correlation data and the mean and variance of the bulk 
depth distribution obtained from linear fuel array transects. The location of cell 
Y¥ nzellative to cells X35 X,5 and Xx. in the cell filling model is: 


Lge 


where the dependent variable, Y, is the cell being filled, and the independent 
variables, X)> X, > and X, are the cells already filled. It is important to dis- 


tinguish between the data ccllected from the linear transects and the requirements 
of the cell filling algorithm. 
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Constants are added to the denominator of each ratio to avoid 
dividing by zero. Appropriate corrections are later made upon 
assignment to the cell. Metric conversion of these histograms 
is not given because the histograms are working tools and do 
not lend thenselves to easy conversion. 


CUMULATIVE PROBABILITY 


0 0. 02 0. 04 0. 06 0.08 0. 10 0. 12 0.14 0. 16 0.18 


10 h LOAD | (BULK DEPTH +1) (Ib /ft?/in ) 


CUMULATIVE PROBABILITY 


1h LOAD | (10h LOAD + 0. 08377 ) 


CUMULATIVE PROBABILITY 


0 0.4 0.8 12 1.6 2.0 2.4 2.8 3.2 


100 h LOAD / (10h LOAD + 0. 08377 ) 


Figure 9.--Cunulative probability distributions for asstgning 1h, 10h, and 100h fuel 
loads to each cell given the bulk depth. 
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The principles of multiple regression analysis are used to estimate Y from 
X,> and X,. 


The following analysis is simplified by transforming X's and Y's to deviations 
from the mean: 


Xi > 


eX eX 
and 

yraryy Shye 
Normal set of equations for Y = £(X, .X, 5X.) + eis 


2 
Dex by t xX Xibo + Xa %sibs = WG 
5%) * yi aiPa + DXoa%ss = LiXi%G 
W%35%145) + 0X34 252 PONS gids = D%3i%i 


and 
DY; = DOG 4b) + Xpgby + Xyyb3) = nby- 


Note that the variables in the fourth equation are uncorrected for deviations from 
the mean. 


The serial correlation coefficient obtained from the bulk depth transect is 
related to the product mean and has this relation: 


a 


COM Geen) 
tia cor(x.,x. | dy = =| 
aS SY Wictel 2s Wane 3 ) 


a 


where a is the lag and b is the orientation and have values a = 1, 2; b = 1, 2, 3. 
The above expression can be written; 


Thus 
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x's and y's are interchangeable (see p. 13) as we move through array filling cells. 
Thus 


and 


a ben ee 
Escl = sins . 


The following diagram is used to determine orientation and lag of the correlation. 


b 


-2 
b=1 GAN) b =3 
RY 
(8) 


As examples, it can be seen that the correlation of cells x, and x, have orienta- 


2 3 
tion 1 and lag 1. The correlation of cells aT and A have orientation 3 and lag 2. 
Cells A and B are included to help in evaluating x,.x The correlation, X14 X35 > 


ila asia Y 
is estimated by averaging the correlations of lag 2 illustrated by cells A and B 


in relation to cell X)> because a direct estimate from the linear transect does not 
exist. 


Since cells A and B are adjacent and directly above and below X 


ee 3? we define 
X14%si as: 


X14%3q = 2] + Tp3)/2- 


Also, we can make the following observation: 


and 


We can now transform the coefficients of the b's of the normal set of equations 
with information gained from serial correlations along the three principle axes. 


The transformed normal set of equations are below: 


— — r +7 — —_ 


2 2 Pil DEN 7? s D 
x Beene and eh (Saper)s Deane ie 
r13x by oP OK b, + T14x be = 1) 5x 
nga) GM — —= —= 
Zi DEN? = 2 
(Geers) Sty Pg 15% 
and 
by = (Ole S b, - b, - bz) 
where 
x = estimate of the mean bulk depth 
ne = estimate of the variance of the bulk depth. 


2 , : ere : : 
Note that although x” cancels in the above equations, it is retained for clarity. 


Transects oriented in directions 1 and 3 should have equivalent parameters be- 
cause of their symmetry about the uphill orientation (orientation 2). The following 
changes have been made to insure this symmetry: 


Fi = Sig ns ele 


Primes pertain to transformed new values. Theoretically, the serial correlation of 
lag 2 should be the square of lag 1. To insure that the two matrix elements (1,3) 
and (3,1) that have correlation of lag 2, (R54 + Ri 3)/2, are influenced by this re- 
lation, the following change has been made: 


et ao ae len 4 et ‘ aval = [ry * 1)3)/2]" 


where the second term is the mean of the squares of the three correlations of lag l. 
Note that Ty and TI3 are the new values as expressed above. 
Solving the three simultaneous equations gives: 


2 


by = b, - R,Q - Tr, )/(R, - 2k, + 1) 
b, = Sib 2R)b) 
where 
Sy Cine at 
5 2 2 
R, = (5) + T,3)/2 + (2R, + r),)/6. 
Rewriting for completeness 
De = MG = = by > DD) 


where X is the mean bulk depth. 


33 


The coefficients of the regression equation, by> b5> and bz, depend only on 


the serial correlation coefficients of the slash transects (all three orientations 
at lag 1 and orientations 1 and 3 at lag 2). The coefficient De has an additional 
dependency on X, the mean bulk depth. 


APPENDIX V: ERROR TERM 


The error term is derived by random access of the cumulative distribution of 
the bulk depth (fig. 10). Random access is achieved by entering on the y-axis of 
the distribution with a random number from 0 to 1. The mean bulk depth is then 
subtracted from the accessed bulk depth to obtain the first estimate of the error, 
e'. The error is further modified by taking into consideration the relation of 
the variance of the estimate of y, y, given x, to the variance of y (Kendall and 
Stuart 1967): 


2 MD 
Sa =) (bh SA ys 
y |x ( sy 


where R = multiple correlation coefficient. 


We can write the following expression for the error by recognizing the simi- 
larity of the variance to the square of the error: 


e = V1 - R? ens 


CUMULATIVE PROBABILITY 
CUMULATIVE FREQUENCY 


0 4 8 12 16 20 24 
BULK DEPTH (in) 


Figure 10.--Cumulative probability distribution of the bulk depth. 
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